CRBLASTER: A Parallel-Processing Computational Framework for Embarrassingly-Parallel Image-Analysis Algorithms

The development of parallel-processing image-analysis codes is generally a challenging task that requires complicated choreography of interprocessor communications. If, however, the image-analysis algorithm is embarrassingly parallel, then the development of a parallel-processing implementation of that algorithm can be a much easier task to accomplish because, by definition, there is little need for communication between the compute processes. I describe the design, implementation, and performance of a parallel-processing image-analysis application, called CRBLASTER, which does cosmic-ray rejection of CCD (charge-coupled device) images using the embarrassingly-parallel L.A.COSMIC algorithm. CRBLASTER is written in C using the high-performance computing industry standard Message Passing Interface (MPI) library. The code has been designed to be used by research scientists who are familiar with C as a parallel-processing computational framework that enables the easy development of parallel-processing image-analysis programs based on embarrassingly-parallel algorithms. The CRBLASTER source code is freely available at the official application website at the National Optical Astronomy Observatory. Removing cosmic rays from a single 800x800 pixel Hubble Space Telescope WFPC2 image takes 44 seconds with the IRAF script lacos_im.cl running on a single core of an Apple Mac Pro computer with two 2.8-GHz quad-core Intel Xeon processors. CRBLASTER is 7.4 times faster processing the same image on a single core on the same machine. Processing the same image with CRBLASTER simultaneously on all 8 cores of the same machine takes 0.875 seconds -- which is a speedup factor of 50.3 times faster than the IRAF script. A detailed analysis is presented of the performance of CRBLASTER using between 1 and 57 processors on a low-power Tilera 700-MHz 64-core TILE64 processor.Comment: 8 pages, 2 figures, 1 table, accepted for publication in PAS

Orbital Phase Variations of the Eccentric Giant Planet HAT-P-2b

We present the first secondary eclipse and phase curve observations for the highly eccentric hot Jupiter HAT-P-2b in the 3.6, 4.5, 5.8, and 8.0 μm bands of the Spitzer Space Telescope. The 3.6 and 4.5 μm data sets span an entire orbital period of HAT-P-2b (P = 5.6334729 d), making them the longest continuous phase curve observations obtained to date and the first full-orbit observations of a planet with an eccentricity exceeding 0.2. We present an improved non-parametric method for removing the intrapixel sensitivity variations in Spitzer data at 3.6 and 4.5 μm that robustly maps position-dependent flux variations. We find that the peak in planetary flux occurs at 4.39 ± 0.28, 5.84 ± 0.39, and 4.68 ± 0.37 hr after periapse passage with corresponding maxima in the planet/star flux ratio of 0.1138% ± 0.0089%, 0.1162% ± 0.0080%, and 0.1888% ± 0.0072% in the 3.6, 4.5, and 8.0 μm bands, respectively. Our measured secondary eclipse depths of 0.0996% ± 0.0072%, 0.1031% ± 0.0061%, 0.071%^(+0.029%)_(-0.013%), and 0.1392% ± 0.0095% in the 3.6, 4.5, 5.8, and 8.0 μm bands, respectively, indicate that the planet cools significantly from its peak temperature before we measure the dayside flux during secondary eclipse. We compare our measured secondary eclipse depths to the predictions from a one-dimensional radiative transfer model, which suggests the possible presence of a transient day side inversion in HAT-P-2b's atmosphere near periapse. We also derive improved estimates for the system parameters, including its mass, radius, and orbital ephemeris. Our simultaneous fit to the transit, secondary eclipse, and radial velocity data allows us to determine the eccentricity (e = 0.50910 ± 0.00048) and argument of periapse (ω = 188°.09 ± 0°.39) of HAT-P-2b's orbit with a greater precision than has been achieved for any other eccentric extrasolar planet. We also find evidence for a long-term linear trend in the radial velocity data. This trend suggests the presence of another substellar companion in the HAT-P-2 system, which could have caused HAT-P-2b to migrate inward to its present-day orbit via the Kozai mechanism

Wiyn open cluster study. VIII. The geometry and stability of the NOAO CCD mosaic imager

The field of direct imaging is experiencing a considerable growth in the number of available CCD mosaic imagers, especially on large telescopes. To fully exploit the astrometric potential of these imaging devices, we develop a technique, utilizing an astrometric standard, for precision transformation of pixel coordinates into a global coordinate system. We have constructed a new astrometric standard set of 1863 stars in the field of NGC 188 to derive the CCD chip constants for the NOAO CCD Mosaic Imager. The multiple-epoch data on the Mosaic's metrics indicate that this CCD mosaic device may have experienced a one-time nonelastic expansion. We also present a new determination of the pixel scale and the optical field angle distortion constants for the KPNO Mayall 4 m telescope prime focus field corrector. To establish a reliable history of the CCD mosaic imager metrics for current and potential future astrometric applications, we recommend obtaining astrometric calibrations for CCD mosaic imagers on a regular basis. Apart from mechanical positioning of the CCD mosaic camera on the telescope, noticeable changes in the thermal environment of CCD mosaic chips should also prompt new astrometric calibrations. It is shown that, following all precautions, the NOAO CCD Mosaic Imager can produce excellent astrometric results on the Mayall 4 m telescope.SCOPUS: ar.jinfo:eu-repo/semantics/publishe